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By Irancis B . . Hildebrand ■ 



SUMMARY 

A least-squares procedure adapted, 'to numerical 
calculation is presented for fhe . approximate solution 
of the Prandtl lifting-line equat ion. Sufficient data 
are tabulated to permit a solution of the equation by 
purely numerical methods for an arbitrary symmetrical 
variation of the chord and the angle of attack." In. 
addit ion,- modified procedures are formulated for the . 
analysis of vings in which the spanwise variation of 
the chord or angle , of attack is discont inuous t ."She 
methods proposed are illustrated by explicit numerical 
analysis of rectangular vings without twijst and with .. 
linear a-nd ■ quadratic twist, a tapered wing with rounded 
tips and partial— span flaps-, and a rectangular wing with 
a cut— out. Comparisons are. made- in several cases with 
the results of other procedures. ■ 

The computation, involved in the procedure, is en- 
tirely mechanical and is conveniently carried out .on. a 
computing machine. She accuracy attained in. a solution 
using the tabulated data should be comparable to the 
accuracy of the given wing data in all practical cases , 
while the time required is considerably less than that 
required bymore elaborate procedures, such as those of 
Lotz and 3etz, and only slightly greater than that re- 
quired by less- exact methods^ such as those .of Glauert 
and Bani. While the modified analys is applicable to a 
wing with a discontinuous angle of attack or chord re- 
quires a small amount of additional computation ,_ it is 
probable that , the resultant accuracy in such cases could 
be attained by the JLot-z procedure only after, a very 
lengthy series of calculations. 
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INTRODUCTION 



The distribution of lift over. the span of a wing, 
in uniform motion is determined, according to the Prartdtl 
theory of the lifting line, as the solution of .a singu- 
lar integro— different ial equation the mathematical com- 
plexities of which are such -that "exact solutions have 
"boon obtained only in very special cases. While several 
methods have "been devised for obtaining approximate solu- 
tions to this equation, it is felt -that a new procedure 
based on a method of least squares which was presented 
in reference 6 may be of practical interest. 

In the usual procedures an approximation to the lift 
function is assumed .as the sum of a finite number of ap- 
propriate approximat ing functions with undetermined coef- 
ficients, after which the coefficients are determined in 
various ways so that the lifting— line equation is approxi- 
mately satisfied. "While it might bo expected that the 
determination of these parameters would be most efficiently 
accomplished by a method of least squares, tho only ap- 
plication of such a method known to the writer (reference 2) 
was not well adapted to numorical computation for arbi- 
trarily varying chord and angle of attack. In addition, 
tho single case treated was that of a wing with discon- 
tinuous angle of attack,, for which tho procedure as^ given 
in referonco 2 failed to give satisfactory results. 

The purpose of tho pr esent - pap or is to present a . v; 
lea's t— squares procedure in which the major part of the' 
numerical calculation can be readily carried out on a 
computing machine, and in which the amount of labor in- 
volved is not dependent upon the nature- 6f~T;he variation 
of the chord and the angle of attack. Since all the 
previous procedures are notably inadequate for the analy- 
sis of -wings with discontinuous spanwise variation of 
angle of attack or chord, an explicit treatment of such 
cases is included. . 

This investigation, conducted at the Massachusetts 
Institute of Technology, was sponsored by, and conducted 
wp.th financial assistance from, the National Advisory 
Committee for Aeronautics'. 
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STMBOLS 



spanwise ooordinate in units of the half span 
(y = 6 at the root, y - ±1 at the tips) 

span 

chord 

root chord ( c( 0) ) 

chord divided by root chord (c/c^) 

angle of attack 
dynamic pressure 

section lift (per unit length along span) 
section lift coefficient ( I f \c) 

auxiliary lift function ( ccj/mc^ , ^/q.™c E ) 
profile constant (dcj/da) 
diEiensionless constant (mc-^/413) 
projected wing area £ ^- f c(y) ±y) 



total lift 




coefficient of total lift (L/q.S) 



angle of attack corresponding to approximate lift 
distribution 
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MATRIX NOTATION 



So presentati o n. of seta of linear equations .— In this 

paper a sot of linear equations is represented in a con- 
densed form. Thus, for example, the formal equation 



A 0 


±1 






a oo 


a 10 






a 01 








a 02 






tt 3 



will "bo used to represent the coefficients of three 
linear equations in the two parameters A 0 and A xl 

of which the first is 

a^ 0 A 0 + a 10 A^ =. a 0 

M atrix multiplication .— A rectangular array of 

elements is called a matrix. 3Tor the present purposes 
it will ho convenient to define two types of matrix 
product a » ■ • - 

Tho first type, which will be called the aJloU 
pr oduc t . is illustrated "by the following example: 



°0 




8*00 


a 10 




c o a oo 


Co a 10 


c l 


* 


a oi 


a n 




c i a oi 


°i a : 1 


C 2 




a 02 


a l2 




C 3 a 02 


c a a i2 



In general, the star product of a one— column matrix^ of 
m elements into a matrix having m rows and n 
columns will be defined as a matrix having m. rows and 
n columns, wherein each element is tho product of tho 
corresponding element in the original m by n matrix by the 
element of the ono— column matrix which lies in its row. 



The definition of tho dot product of two matr-ices 
may bo illustrated by the following example: 
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b oo ;Ko 



J 01 
3 0B 



J ll 

3 is 



<2T 



?oo. 

a 01 
^03 



''lO 
'11 
L 12 



(^oof-oo-^pif^^s^oe)' ^oo^io+^oi^i + ^ceais ) 
C^io aoo+^x a 0 i + ' h is a oB^ &io *h.o +h ii a n +l3 is a m) 



In general, the dot product of. two rectangular matrices 
having an equal number:- of rows will "be defined as a 
matrix wherein the element in the ^—th row and the k— 4h~ 
column is the algebraic sum of all pr oduct s _ of ^ corr e^- 
spending elements in the J— th column of the first factor 
matrix and the k— th column of the second factor matrix. 
(This product is equivalent to the conventional matrix 
product of the transpos-e of the first matrix into the 
second matrix.) 



SOLUTION OE-THE L IET ING— L INE EQUATION M A METHOD 

. . \. -* - ' — «■ 

': OS 1 LEAST SQUARES 



According to the lifting— line theory, the lift 
l(y) per unit span acting on an airfoil is determined 
in form's of the chord- c(y) " and -the angl e., -of, attack ' 
a(y~) '-by the intogro— differential equation : 

Uy) = c*ly){i*<i«£?l^ vJ^W^ Ua> 

and tho "boundary conditions -" .. 

I(±l) =0 , ,.. ... . (lb) 

whoro y is a spanwise coordinate measured from the 
root in units of tho half span Jb/2, c*(y) is the 
ratio of the chord to the root chord, 

'-: c *(y) = . (<>?) 

. ■ • c 2 • . .. - 

4 ■. 'is .-thotdynamic pressure , m is a .profile constant, . , 

dot \ '• ■ • ... .•f. ■ «, " 

- i , ' and u 'is" the dimensionless constant defined 

day - -. • T-- 

by the equation " . ' ' 
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Tho integral' appear irig.. in equ.atio.n- (l) is not a proper 
intogral, but is to "be assigned its Cauchy principal 
value, according to the definition .. . 



f- ^o{f_ x + f + J 



If the notation 



3P(y) = JUL ' ; (4) 



is introduced, equation (l) can "be written in the form 



c*(y) 



dF dn 

dTl y-Tl 



- a(y) 



r(±i) - o 



(5) 



Th-e usual method of obtaining . an approximate solu- 
tion to this equation consists in assuming that F(y) 
can "be approximated satisfactorily in the range [yj^ 1 
by a finite series of appropriate functions 



F(y) « ) A a cp fe (y) 
n^ 



so that equation (5) becomes 

H 



where 



n = 1 

•n(y) - + trj^-nr y^T 



:(6) 



(?) 



(8) 



and in then determining' the parameters A n by requiring 

that equation (7) b© a true equality at H arbitrarily 
chosen points (reference l). This procedure thus leads 
to a set of N linear equations 



r 



in the 
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J A S ®n kip =a(y k ) 
k» l, 8, ... I 



undetermined parameters A a . 



(9) 



A more elaborate procedure consists in determining 
the parameters "by requiring that the integral of the 
square of the difference "between the tws sides of equa- 
tion (V), over the range . [ y J ^ 1, he a minimum .( ref- 
erence 2) , •' . * 

_ N _ 2 



n = l 



. dy = minimum 



(10) 



Equating the partial derivatives to zero leads again to 
a set of IT linear equations 



) &nf «*(y)*a(y)47 = £ * p (y)&(y)dy^ 

. . . B" J 



n= 1 



(11) 



p = 1, 2, 



which serves to determine the constants '. A n . 



squat 
equat 
here 
cedur 
other 
basic 
by a 
equat 



In order to avoid the integrat 
ions (ll) and reduce the appro 
ion (5) to a purely numerical 
presented a. modif icat ion of th 
e which has been used previous 
problems! (reference 6). This 
ally in approximating the inte 
sum of weighted values of the 
ion (10) is replaced by the co 



ions involved in 
pcimate solution of 
process, there is 
e least square pro— 
ly in connection with 

method consists 
gral -in equation (10) 
integrand, so that 
ndit ion 



fir 



& 



= minimum 



(12) 



^ Dk [ 2L A a'*a(3 r k)- a (yk)] 
k^i n=l ' ™ • • 

where Djj. is the integration coefficient associated 
with the value of the integrand at the point y^ . . 

Equations (ll) are now replaced by the following set of 
N linear "equations: 
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N M .% .. 

7 A n ^ E k $ p (y k ) $n (y k ) = ^ D k $ p ( y k ) a ( y k ) 



n = l k=l 



k = l 

p = 1, 2, . . N 



(13) 



With the notation of the preceding section, the 
coefficients of this set of epilations can be written in 
matrix form 

Dk*a(*k)l e|« a (7k)| = ||Dkt a (yk)l"o||a(irfc5|| ("> 

where the row and column indices are k and n, re- 
spectively. It follows that the coefficients in equations 
(13) can he conveniently obtained by first writing the 
set of equations ■■ i . . ' 

N 



y An *n(yk) = a(y k ) 



n=l 



k - 1, 2, 



M 



(15) 



in matrix form, and by then forming the dot product of 
the auxiliary matrix 

into "both, sides of the resultant formal equation. 

While equations (15) are of the same form as equa- 
tions (9), the points y- k now correspond to the weighted 

ordinates in equation ( 12) , so that the range of k is 
arbitrary. Thus, in place of satisfying equation (7) . 
at a number of points, N, equal . t.o the number of un- 
determined parameters, the present procedure satisfies 
this equation as nearly as possible at an arbitrary 
number of points, M . 

Since equations (13) are homogeneous in. the weight- 
ing coefficients D k , these coefficients may be chosen 
as any convenient multiple of the actual integration 
coefficients. Also, it may be desirable in certain 
cases to interpret equation (12) otherwise than as an 



F 
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approximation to equation (10). That, i.s , in place of 
requiring that the mean square -of the difference "between 
the two sides of the lifting— line equation "be a minimum, 
it may be desirable to require 'that the equation be more 
nearly satisfied at some points than at others. The 
coefficients Djj. are then of the nature of "influence 
coefficients" and may be determined by magnifying cer- 
tain of the integration coefficients in proportion to 
the dogree of satisfaction desired at- corresponding in- 
tervals along^ Vhe : _span. . 

If .the computed values of the parameters are sub- 
stituted in the" left— hand sides of equat ions- ( 15 ) the 
values of the angle of attack 5 corresponding to the 
approximate lift distribution specified by equation (6) 
are determined at • M points along .the span.. A compari- 
son of this angle— of— attack distribution with the pre- 
scribed distribution a will give an indication of thje 
degree of approximation attained in the solution. 

, OUTLINE OP THE PROCEDURE POR NUMERICAL SOLUTION 

OP THE EQUATION, . . 

In this section a procedure involving a five— term 
approximation to P(y) ■ is ■ expl icitly developed for the 
case when ' c(y.) ■ and a(y) are symmetrical with respect 
to tho wing root. An analogous procedure can be developed 
for the treatment of the case of an ant i— symmetrical 
angle of attack. • 

It is convenient to choose the approximating fuEc- 
tiohs. tp n of equation (6) so that the functions defined 

by tho integrals- • -• 

iL r 1 dcp s ^ 

*T -J_i dn y-ri 

have a simple analytical expression and are regular in 
the neighborhood of the wing tips (y = ±l). Such func- 
tions are, in general, characterized by the property 
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BP- -I 

9 n (y) ~ (const-. )' (l-y s ') 2 

y -r*±-l 



(16) 



where n is a positive integer. In addition, .the 
functions .tp^ should he of a form readily adaptahle 

} 

to the approximation of the function 3 1 , the character- 
istic behavior of which usually is known. 

In the present procedure an approximation to 3T(y) 
is assumed in the form 

3 

T(y) *B : y.«log +vT="y^ V A an 7» a (?•?> 

n = <5 

->--w.With the exception of the first term, the approxi- 
mating functions are conventional ones employed else- 
where, Ihe coefficient of B, which is of the form 
required by equation (16), was originally chosen for use 
in cases when a(y) has a discontinuous first dorivativo 
at tho root (e.g., in the case of a symmetrically linear 
angle of attack), since the contribution of this term to 
the integral representing the induced angle of attack, 



has a discontinuous derivative at tho root (y=o), whilo 
the function itself has a continuous derivative at this 
point. The function was, however, retained for use in 
the mora general case since it complements the other 
approximating functions, being intermediate in behavior 

Wfeweefi the functions </l— y» and ys </l- y s . (-See 
figure 1, ) ■ 

With the approximation of equation (17) the lifting- 
line equation (5) becomes 
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r ft.a(y) [r\ 1 (18) 

•3?or the approximate integration" -indicated in equation 
(ll) nine points y k are chosen, equally' spaced over 

'the interval O^y-^l*; so that ' y k = 1z.fi , 6 : 

*k = 0,1, '2, . . . 8. -'If equatie-n (18) -is evaluated at 
these points, nine linear equations in the five parame— 

'ters . ar e. .obtained the coefficients ef which. .are written 
in matrix form in equation ( 19 ) . '!■'•.< 
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1 

0*(i/4> 

c*U/fc) 
o*(i/s) 

C*U/ 8 ) 
C*( 3/4) 



11 



B 

.00000 
.04326 
.12897 
.23018 
.329 24 
.40897 
.44739 
.40452 



A 0 
1. 00000 
.398216 
.96825 
.92702 
.86602 
.78062 
.66144 
.48412 



Lc*( y )J y=1 Lo*(y)J ys=1 |_c* 



.3 
. 00000 

.01550 

.06052 

.13036 

.21651 

.30493 

.37206 

.37066 

Vl-y 3 " 



ty). 



A* 

.00000 
.00024 
.00378 
.01833 
.05413 
.11911 
.20928 
.2837.8 



A* 

.00000 
.00000 
.00024 
.00258 
.01353 
. 04653 
.11772 j 
.21727 



7 = 1 



[^EzS BEES: 



B 


A 0 


A 2 


A 4 


A* 






-1.00000 


1.00000 


-.50000 


-.12500 


-.06250 




a( 


o) 


-,60730 


1.00000 


-.45312 


-.14722 


-.06894 




a( 


i/e) 


-.21460 


1.00000 


-.31250 


-.19922 


-.09309 




a( 


1/4) 


.17810 


1.00000 


-0 07812 


-.23706 


-.14521 




a( 


3/8) 


J. .57080 


1.00000 


-J85000 


-.18750 


-.20312 




co( 


*/■) 


.96350 


1.00000 


.67188 


.05200 


-.17322 




a( 


5/8) 


1.35619 


1.00000 


1118750 


.61328 


.181S0 




a( 


3/4) 


1.74889 


1.00000 


1.79688 


1.65747 


1.34811 




a( 


7/8) 


2.14159 


1.00000 


2.50000 


3.37500 


4.06250 




a( 


o 



(19) 
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The elements in the Last row *f-' the second matrix are 
all zeros, unless the chor d ts psrs to zero at the tips 
at least "as rapidly as. vl-ys ,: -In case the wing ap- 
proximates . an elliptical plan f-orm near the tips, all 
thes e ' 61 ement s ,have an eq.ua.1 non— zero value which can 
be determined in a : simple manner. If a /constant X is 
determined so that , in the immediate neighborho od' of the 
tips, ... 



'720) 
1/X. 



then these elements are to Tfe assigned the value 
While cases in whi'ch the^ chord tapers to zero more 
rapidly than ■ -Jl — ys ■ 'involve certain mathematical dif- 
ficulties, it is: probable 'that for wings with rounded 
tips the' error. ^ introduced by us ing a tip approximation 
of the form given in equation (20) would not be great. 

If now the values of the chord and angle of attack 
are known at the nine points^ along the span, and if n 
is prescribed, equation '( 19") 'can be put_ in the form 



B 






T>o 


■ a'oo 


" :;a so 


*i 


a oi 


a. Bl 


■ 




* 


* 




- < * 






■ 


^8 


&os 


as e 







a o 






a i 


* 


St 




■* 
* 




CX 8 



(21) 



where, for example, 
b 0 = -[i 



b, = , /„ ^ — .60730^ 



.06052 n. nr.* 



^T787 r . . 



+ 4.06250 |i 



a k = a- (k/8 ) 



(23) 
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In accordance with the results of the preceding 
section, the coefficients of the linear equations (13) 
are obtained "by forming the dot product of an auxiliary 
matrix' into both sides of formal equation (2l). This 
auxiliary matrix is formed by multiplying the elements 
of the k — th row of the left— hand matrix of equation 
(21) by the weighting coefficient associated with 

point y^ in equation (12), If it is required that the 

mean square of the difference between the two sideB of 
the lifting— line equation be as small as possible, so 
that equation (12) is to approximate equation (10;, 
these coefficients are to be proportional to a Vet of 
integration coefficients, The best results have been 
obtained, in' general, if coefficients proportional to 
those of Simpson's rule are used (for a detailed dis- 
cussion, see reference 6, pp. 319— 523), so that 

Do - D 8 * 1/2 
D 1 « D 3 = P 5 = D 7 = 2 „ 

D 2 = D« - » e = 1_ 

and the auxiliary matrix is »f the form 



(26) 



If the dot product of matrix (25) Into both Bides 
of equation (21) ±i formed, the coefficients of a set of 
five linear equations in the five parameters B, A 0 , . . . . A 6 

are obtained the solution of which, in connection with 



(24) 

> 



*>0 


£ a 0O 






Sb^ 


2 a 0 i 














» 

■ •* 




• 


V 

• 


4*8 


* 


¥ 





r 



t 
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equation (l?), determines the lift distrTbut ion over the 
span of the wing. ' Ihe amount of numerical work involved 
is considerably reduced if use is made of the fact that, 
as can "be shown,' the matrix, of coefficients in the final 
set of equations is symmetrical with; respect to its 
principal diagonal-. -_ - — 

The data required to compute t^ie values of F(y) 
at nine points are tabulated in' the second matrix of 
equation ( 19 ) (if the elements of the last row are re- 
placed "by. zeros.). Thus 4 f or ■ examp 1 e , . 



1(0) = A ( 



P(l/8) =* 0,04326B * 0.99216A 0 + 0. 01550 A a + 0, 00024 A 4 •> (26) 



= 0 

With the auxiliary function I*(y) known, "the lift dis- 
tribution l(y) is determined from the equation 

"' ■ l(y)' = m .q; c R 3P(yj 

and the section lift coefficient cj(y) follows fr'um 
the equation . 

c(y) 

If L is the total lift acting on -the wing and S 
is the projected area of the wing, the coefficient of 
total lift, Ox,, defined "by the relationship 




is given- by the formula 

°L = -i-y 3P(y) dy 
or, with the approximation of equation (17), 
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°L - ^-y-[i B + A ° A a+ .f.A 4 + :fj A 6 ] (27) 

An indication, as to the accuracy of the solution 
is obtained if the .values of the angle of- attack a 
given by the left— hand sides of the equations repre- 
sented "by formal equation (2l) are compared with the 
prescribed values of a; for example, 



a k = *k B + a ok A o + ' • • + B 6k A s~a k (38) 

It may he mentioned that in place of using the ap- 
proximation of equation (17), the lifting— line equation 
may he first transformed by the substitution 

y = — cos 9 o^e^n 

after which, if an approximation to the function J?(y) 
is assumed of the form 

S , (y)- w B(cose 9 )logf ) ,+ VAn sin n9 (17a) 

\ 4' cos 0j / L-> 

equation ( 18 ) is replaced by 

+ n ^ B(n fcos ei — L) + Y n A n ~f^g^| w a (cos 6) (18a) 

and the present least— squares procedure is .again applicable. 

In the preceding developments it was assumed that 
the profile constant m does not vary along the wing 
span; that is, that parallel chord sections of the wing 
are geometrically similar. If this is not the case, the 
effect of varying M can be taken into account if in 
equation (5) and the' subsequent equations the function 
c*(y) is replaced by the function 

G *(y) 

m. a 



NACA le-chnie-al Note No..' 925 



17 



and the definitions of. the constant, ' and the funt 
tion £(y) are modified -so as to read' 



b 



(3a) 



E(y) = 



V(y> 



4 m R c R 

where m^ is the value of ia at the root. 



(4a) 



EXPLICIT SOLUTION EOR A RECTANGULAR VINC 



In order to illustrate the procedure, of the pre- 
ceding section, a wing of rectangular plan form and as— 

"b 

pect ratio- — = m (fe6) is consfdered. In this case, 
equations (.2) and (3) "become*. ' ■ •• 



and 



>(y>- = i 



(29) 
(30) 



Each element of the left— hand matrix of equation 
,(2l) is thus the sum of the corresponding element of the 
second matrix of equation (19) and one— fourth the corre- — 
spnacLing element. of the third matrix of that equation, so 
that equation (21) "becomes . . . . _ .... 





' A 0 


A 3 


A 4 


A 6 




.25000 


1.25000 


-.12500 


-.03125 


-.01562 




cc 0 


.10856 


1,24216 ' 


-.09778 


-,,03656 


01723 




a i 


.07531 


1.218 25 


-.01761. 


-.04602 


-.Q2304 




a s 


.27470 " 


1.17702 


-.11083 


-.04093 


-.03372 




00 3 


.47194 ' 


.1 .11602' 


.27901 


.00725 


-.03725 




CG 4 


.649 85 


1.03062_ 


.. . .4729 0 


" ,13211 


.00322 




a 5 


.78644 


.91144 


.66893 


,36260 


.16307 




a 6 


,84175 


.73413 


,81988 


.69815 


,55430 




a 7 


.53540 


.25000 


.62500 


.84375 


1.01562 




0l>8 



(31) 



1 
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The auxiliary matrix -is "now formed by multiplying 
the elements in the k— th row of the left-hand matrix 
by the weighting coefficient. Dfc, as given in equation 

(24). If the dot product of this matrix into both sides 
of equation (3l) is taken, the following five linear 
equations are obtained: 



B 

3,47075 
4.19825 
2. §1636 
1^84739 
1.30533 

where 



Ao 

4.19825 
13.43179 
3.09.615 



2.91636 
3.09615 
2. 5 6410 



A 4 
1 .84739 
1.47864 
1.77884 



1.30533 
.89447 
1.32531 



1.47864 1.77884 1.50586 1.26747 
. ,89447 1.32531 1. 26747 1.16176 



a 3 
a 5 



(32) 





-.'12500 


,68500 


-.06250 


-;01562 


-. 00781 








-.21713 


2.48431 


-.19556 


-.07312 


-.'03446 




a i 






. .07531 


1,21825 


-.01761 


-.04602 


- : . 02304 




a s 


a 3 




.54940 


2.35405 


.22166 


-.08189 


-.06745 




a 3 


&z 




.47194 


1,11602 


.279 01 


,00725 


.-,03725 


8 


0-4 


a 4 




1.29999: 


'2.06125 


.94580 


.26423 


, 00645 




ttfl 






.78644 


.91144 


,66893 


.36260 


.16307 




cc.e 




1. 68 349 


1.46825 


1.63975 


1.39630 


1.10860 








.26770 


.12500 


.31250 


.42188 


.50781 




C&8 



(33) 



It may be remarked that the calculation of each 
element in equation (32) involves only a single continu- 
ous operation oh a computing machine. S*or example, the 
coefficient of A 3 in the second row of equation (32) 

was obtained as the algebraic sum of products of corre- 
sponding elements of the second column of the auxiliary 
matrix (the first matrix on the right— hand side of 
equation (33)) and the *hird column of the left— hand 
matrix of equation (3l).: 



(1.21825)( .01761)+.... 



-( .62500 )( .12500) - ( 2 .48431 )(. 09778 ) 
. . . + ( .12500)('. 62500) ='3.09615 

Since the matrix on the left— hand side of equation (32) 
is symmetrical, it was not 'necessary t o compute the 
elements below its principal diagonal. 
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The solution of equations (.32) is also conveniently 
carried out on a computing machine., the process "being 
considerably shortened if account is taken of the symme- 
try in the matrix of coefficients. (See, for example, 
ref er enco 5 i ) 

In the case of a ' uniform angle 1 of attack , for which 
a(y) = aj = constant, the right—hand matrix of equation 
(32) as determined from equation (33), is found to be 

1 4.79185. . 
12.3635.7 
3.79199 
'2.23562 
1.61592 

and the solution of equations (32) then gives 
B = —0.25469 arjj A 0 = 0.82881 arjj A 3 = 0.98119 

A 4 = 32373 A 6 = 1.36382 arp 

The coefficient «f total lift 0 L is determined "by 
equation (27), 

e L = 0.73065 M arji (35) 

In table 1 and figure 2(a) the auxiliary lift func- 
tion #(y) (computed from equations (26)) is" compared 
with the corresponding solution given, by Glauert (refer- 
ence l). This solution was obtained. by assuming a four- 
term approximation to !(y) of the form given in equati 
(17a) (omitting the first term) and by then satisfying 
the lif t ing— line, equation exactly at four points, not in- 
cluding the wing tip. The value's of the left— hand side 
of the basic equation (5) (computed from equations (19)) 
are compared with .the corresponding values determined 
from the G-lauert solution in table 2. It is seen that 
the G-lauert solution satisfies the lif ting . line equation 
extremely well except ■ in the . immediat e vicinity of the 
wing tip where a large, deviation occurs, while for the 
pres/enj; solution the equation is reasonably well satis- 
fied along the span. The two lift distributions agree 
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very closely, however , except near the tip of the wing, 
and the G-lauert lift- coefficient, Cj, = 0.729 ta atp , is 

nearly identical withthe result of equation (35). 

In the case of a symffl,etr.ic.al3,y linear. .angle of attach, 
for. which a(y) = arjjyj , -^the right— hand matrix of equa- 
tion (32) is found to be • 



3.57658 
5.43746 
3.03389 
2.033-02 
1.550i8 



and equations (32) give 



B = 0.11229 aj 



A 0 = 0.16731 arp 



A 4 = -1.67588 aj A 6 



= 1.48110 a 



The coefficient of total lift is then 

0 L = 0.33270 m aj 



2 = 1.2505 7 

Hi 



J 



The variation of the function 3T(y) is presented 
in table & and figure 2(b) in cCbmparison with the 
G-lauert solution, while in table 2 the left—hand side of 
equation (5) is evaluated for the two solutions and com— 
pared with the prescribed right— hand side. In this case 
the actual angle of— attack a corresponding to the 
G-Lauert solution deviates appreciably from the prescribod 
angle of attack a both near the root (y = 0) and near 
the tip (y = l).. The root ' deviat ion , which is due to 
the fact, that the curve representing the function a»(y) 
has- a discontinuous derivative at . y = 0, is decreased 
in the present solution by the presence of the first term 
in the approximation of equation. ('17), and the tip devia- 
tion for the present, solution is alee less pronouncod. 
Corresponding differences occur between the two lift dis- 
tributions, and the G-lauert lift coefficient, Cj, = 0.320 m. a^, 

differs from the result of equation (37) by about 4 percont. 
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In the case of a quadratic angle , of attack t for 
which a(y) = aj y s , the~r tght— hand matrix of equation 

(32) takes, the form . 

: i{2.703-25 
..... 3.29199 
am 2. 41 043 j 

-Li 78 4381 

Jj..43b65l 

The solution of equation (32) is found -to he 

B = -0,20524 a T A 0 = b. 06665 arji A 2 = 1.20613 a^. 

A 4 = —1.00759 atjj A s . = 1.13356 aj 

and equation (27) gives 

Oj, = 0.20607 aoj 

The 'lift distribution corr esponding r 't b this solu- 
tion is presented in table 1 and figure 2(c) while in 
table 2 the functions a and a are compared. , This 
problem was not considered by G-lauert . The' lift' dis- 
tribution in this case j.s more nearly concentrated at 
the wing tip and the maximum'.. lift value is only about 
42 percent of the maximum value for the corresponding 
wing without twist. Also, a higher degree of approxi- 
mation is indicated by the agreement of ..the, two . s ides 
of the lifting— line equation. .-■ ,. 

MODIFIED PROCEDURE FOR VIKGS VITE D IS CONT IHTJOUS 

AITG-LE OF ATT ACE 



It is known that at a point of discontinuity in the 
angle of attack the lift— distribut ion curve has an in- 
finite derivative, so that an approximation of the type 
used in equation (17) could not be expected to give an 
accurate result. In such a case, in order that the right- 
and left— hand sides of equation ("5) have -the same dis- 
continuity at a point, the function I*(y) must be con- 
tinuous but must have the property that the integral 
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JL-f 1 J* 
tt J_ L dn 



has the same discontinuity as the angle of attack, & . 

As can "be ver if ied ■ in- an elementary manner, the 
function r ... ... 



\ . Ty-al / 



satisfies the equation 



" 3 J«i dn. 



dn " 
y-Ti 



COS 4 % 

■ TT 



1 — — ■ cos 1 a 

TT 



a< y ^ 1 



and has the crrrect behavior at the tips, since 

r a(y) "* Wl-a 3 </l~y3 y .— > ±1 

Thus the Integral ». . 



C4i) 



is piecewise constant in the interval fy t ^ 1 and has 
a jump of magnitude ttm> . at the. point y= a. Suppose 
that a(y) is continuous except for finite jumps of 
magnitude J n at the points y = a n . 



a ( a n + ) — a ( an—) = J n 



Then, if the function 



(43) 



is defined, it follows that the expression 



r 
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'/ * ■■ 1'-, L. . v ^- • 



is continuous. 



Ehis fact suggests the following jrocedure for 
solving, the lif ting^lipe .equat "ion ■ (;5) -in such p^ses\ If., 
'an ■approximation'' 'to the"' lif f d'is tribut'SoV fiinct ion .is 
assumed of the form 

where • P is-, defined bosquet i« A (43) and -P 0i .is given 

"by the right-hand side of equation then equation 

(5) can .be written ■■- ^ ..... . .. > .i- ... ■: t-.- S -t : " 

P 0 (y) ^l^ fa-iLv;^).: (46) 



c*(y) 17 J_ l dT i J^-n 
where ■ 

*^ ^ n T 1 dP dr, „ ^ , ' T .: : 

a*(y) = — + — 4 — - a,.(y) 145) 

c*(y) " J_, dn y-T) 

Equation ( 46 )" is " formally equivalent to equation (5) if 
a(y) is replaced. J>y a*(y) -and, since _ct*(y) is con- 
tinuous the solution can "be carried out as outlined in 
the preceding sections. . ... J. 

In calculating the total lift coefficient, 0 L r the 
value of the integral 

- ;. - ■ • , ; '. >;> •; 

r 1 r & (y)dy= -''^yi-a 2 ( 49 ) 

is needed. 



'As an example, suppose that the' aileron deflection 
of a wing is such that 

0 -1 s£ y < -a 

a(y) = -l 1 -a « y < a (49) 

0 a < y ^ 1 



r 
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In this case equation (43) becomes 



or 



, . -if/ % Hl-ay + Vl- a s vft—y 3 1 . f n-ayWl-a 8 ' 
P<7) = |;{-(y-a) !o g [^-^Zlj J +(y+a)l °4 |yTa| 



Also, using equat ion ( 40) t it is found that 



i>«/ 1-y 2 



(50) 



n dn. y-n n 8 I J-i dn y-n J_ x dri y-rU 



-£- cos + C~ ~ 



(tt — cos 



-1 



.)!} 



~1« y< - a 



— £-•=■ cos _1 a] + [l - ^- (tt - cos" 1 a)3 ~a< y < a 

* TT TT 

_ [i_JL costal + CX - — (n-cos~ 1 a)3 a <y ^ X 



dr| ^ 

y-n 



2 _, 

y(- cos. * a — 1 

— cos A a— 1 

TT 



-i-l £ y < — a 
-r& *i y < a 
a < y £ 1 



= a(y) - (l - -f- cos 1 a^ 

* ■ 

Ihus, finally, equation (47) gives 

a*(y)=|^-,(l.- |- cos-a) 



(51) 



(52) 



T 
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EXPLICIT SOLUTION -.JOB, A TAPERED WING WITH SOUNDED TIPS 

AND PAST IAL— SPAN P1APS 



As an illustration of the modified procedure, a 
wing analyzed in reference 4 is now considered. The 
variation of the ch-Ord ratio- c(y)/c £ is represented 

in figure 3, while the following additional numerical 
data are. given: ■• - 



H = 0.1919 



= 10 



(53) 



In the immediate neigh "borhood of the wing tips, it is 
found' that, approximately, . ■ . 



c*(y) s 1.64 «/l - y 3 



±1 



and it follows from equation (20) that 



X = 1..64 



(54) 



An' angle of attack of ■ 1 radian' from 0 to 0. 489' and 0 from 
0.489. to- 1 is prescribed, so. that, in the notation of the 
example of the preceding section, .'*".-'" 



a = 0.489 



(55) 



The equation to "be solved is then 



FoCy) + U 
c*(y) 



dP, 



dn 



dTi y— ri 



co*(y) 



Where 



(56) 



Po(y) = P(y) - P(y) 



(57) 



and P(y) and a*(y) are defined in equations (50) and 
(52), 
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The functions c*(y), P(y) and 
at the nine points considered in 



a*(y) are 
table 3. 



tabulatod 



If the operations indicated in equation ( 19 ) (where 
a is replaced by a*) are carried out, equation (2l) 
takes the form 



B 


A 0 


A 2 


A4 


A 6 




-.1919b 


1.19190 


-. 09595 


-.02399 


-.01199 




1.85269 


-.07040 


1.25021 


-.07042 


-. 02800 


-v01323 




1.94784 


.10621 


1.29847 


. 00920 


-.03391 


-.01759 




1.94088 


.31748 


1.33285 


.14545 


-.02293 


-.02469 




1.76657 


.54852 


1.34659 


. 33665 


-.03619 


-.03094 




1.09D23 


.77976 


1.32734 


.57247 


.18323 


,03444 




. 4*026 


,97608 


1.25020 


.82318 


.45254 


.22316 




. 05273 


1.05476 


1,05256 


1.00377 


,82257 


.64496 




-. 17629 


1.02073 


.80166 


1.08951 


1.25742 


1.38935 




-.32527 



•(58) 



The coefficients of the final set of linear equations, 
obtained by taking the dot product of the auxiliary matrix 
defined in equation (25) into both sides of equation (58), 
are then found, 



B. ■ A 0 

5.45683 7.35230 
7.35230 18.51243' 



4.66690 
3.11236 
2.31514 



5.71837 
3.14707 
2.12795 



4.66690 
5.71837 
4.11195 
2.92894 
2.26282 



A 4 

3.11236 
3.14707 
2*92894 
2.42110 
2.05002 



2.31514 
2.12795 
2. 26282 





1.76261 




15.54154 




.62091 




-.51871 




-.607 23 



(59) 



and the solution of the corresponding set of equations is 



B = 4.20498 
A4 



A.o = 1.40506 
6.24691 A 6 



-9 .9 4467 



-1.96339 



(60) 



In table 4 and figure 4 the auxiliary lift function 
I'(y), determined by the equation 

y(y) = P(y) - 3*0(7) 

is presented in comparison with a solution given by Pearson 
(reference 4). This solution was obtained by using a ten- 
term series of the type given in equation (l7a) (omitting 
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the first term) and determining the parameters "by a 
method given by Miss Lotz (reference 3). For further 
comparison there is included in. figure 4 a solution ob- 
tained by the present least— squares procedure with an 
approximation of the type of equation (17), that is, 
without, using the additional approximating function 
P(y). While the Pearson solution agrees closely with 
the present solution over a large part of the span, it 
appears that, even with a t en— t erm ■■ appr oximat ion , the 
correct behavior of the lift curve cannot be satisfacto- 
rily approximated near the end of the flap except by the 
use of a term similar to the function P. The values 
of the angle of attack corresponding to the present solu- 
tion and to the Pearson solution are compared with the 
prescribed values of a in table 5 and figure 5, - 

MODIFIED PEOOEDUEB POfi VIKTGS WITH DISCONTINUOUS 

CHOED YAEIAEION 



Suppose that the angle of attack co(y) is continuous 
and that the chord c(_y) is continuous except for finit.e 
jumps at the points y '= a n . Then the function P~(y)'/c*(y) 
has corresponding discontinuities of magnitude — Y n P(a n ) 
where ■ * 

o*( sn+ )~ °*( a-n-) :.-"(6l) 
c*(a n +)o*(a n -) 



^n=- 



- c*(a n +) 



in")]" 



Hence, since a is continuous, the left— hand side of 
equation (5) must be continuous and the function 

&_ dJ dn . • . 

w J_ x cLn y-n 

must have discontinuities of magnitude +*V n P(a n ) ' at the 
points y = a a . 



If the function F(y) is written in the form 



P(y) = v x M - Q(y) 



(62) 
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whore • 

1 



Q(y) = -7- > Kn r an (y) (63) 

n 

and T a (y) is defined by equation (39), the constants 
K n can he determined so that the__ 1 eft— hand side of equa- 
tion (5). is continuous and consequently the function 
I 1 1 (y) has a finite derivative inside the interval jyj<l. 

Tot, according to equation (40), if equation (62) is in- 
troduced into equation (5) the left— hand side of the re- 
sulting equation has, at each poin^ . y = am , a discontinuity 

of magnitude. 

n - 

It follows that the discontinuit ies" will disappear if the 
constants E n satisfy the equations 



5 K n \ ^m r a* ( «o ) -~ n Sjm [= 17 |i Y m P ( a m ) 

J (64) 
m = 1 , 2 , 3 , ... 



where _ - 

o, m y n 
Sffln^ . • (65) 

I 1 , m = n 
L 

The constants K n are thus determined, from equation 
(64) as linear combinations of-the values of S' 1 (y) at 
the points of: dis coixt-inuity . ■ 

Equation (5) can now he written in the form 

?i(y) 



c*(y) 



+ 



iraii-i%> t ±rMi] =a(r , (66) 

17 «-Li dri y-t) . i_c*(y) n J_i dri y-ri J 



r 



• NACA • 0? e ohh'tis&V Sr» t e ' ifo . ~9"SB 



Since the function ^.(y) has a finite derivative inside 
the interval' jyl <1, an . appr oximat ion to ^(y) can be 

assumed in the' firm given "by equation (17) and the term 
in brackets in equation (66) becomes a known linear combi- 
nation of the parameters B and I f this term is 

evaluated at the points with the help of equation 

(40) and if the- set .of coefficients if B and A 



2a 



in 



the resultant expressions is written as a matrix, the 
matrix equation replacing equation (2l) is obtained by 
subtracting this matrix from- the left— hand side of equa- 
tion (19). The least— squares- procedure can then be ap-^ - 
plied as before. 

As an example, suppose that the chord variation of a 
symmetrical wing has discontinuities of equal magnitude 
and opposite sign at the points y = ±a and write 



V = 



'(a+) - e*(a-) = _ c*(-a+) - c* (-a-) 

c*(~a+) c*(~a-) 



c*(a+) c*(a-) 



(67) 



Then equation (620 becomes 

C(y) = -^j- [z-aF-a(y) + K a r a (y)J 



(68) 



Also, since 3P(y) . is an even function of y and since 



r a (a) = r_ a (- a ) = 0 

r a (-a) = -T_ a (a) = 2a log 



a 



(69) 



equation (64) becomes 

-TTuK-a -*,;^a log a K a = - n (i Y ^(a) 
— 2 V a log a K_ a - tt fj. K a = P M> Y ( a ) 



^ (70) 



so that 
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K a = ~ K -a 



TT (J. 



p x (a) 



(71) 



2 Ya log a — tt|j, ' • 
Thus equation (68) becomes 



Q(y) = — — — y_ a (y) - r a (y)^ » *(y) ( 

' 2 Ya log a-TT(J, ** tt u, - 2 Ya log a 

where P(y) is defined by equation (50), and the in- 
tegral equation (66) oan be written in •fch.e form 



c*( 
wher e 

P(y) - 



— + _ J p( y ) p l ( a ) - a(y) 

y) n «i 1 d 



(73) 



dn y-ri 



TT n Y 



tt(Ji — 2 Ya log a v -c*(y) 



^ P(y) M. p dp dri *\ 
^c*(y) tt dn y — n*' 



(74) 



EXPLICIT SOLUTION POP A RECTANGULAR WING- WITH A CUT-OUT 



As an illustration of the procedure of the preceding 
section, the effect of a rectangular cut— out on the lift 
distribution over a rectangular wing with an original 

b 

aspect ratio of .magnitude — = m without cut— out is 

c . — 

analyzed. The length of" the cut— out is taken to be one- 
quarter of the span and the width to be one— quarter of 
the chord (fig. 6). In this case there follows 



16 



a = — 



e*(y) = 
and equation (67) gives 




Lyl < % 



(75) 
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Y = 



(76) 



Iron equation (50.), 



i + -J-y +^ Vl-y 2 



] 



f* 1 



and from equation (5l), 



- < 

y-Ti 




i_ 

4 




(77) 



(78) 



If ^(y) is assumed in. the form given "by equation (17), 
there follows _.■*■" 

P z (a) = 3* j (.-J-) = 0. 12879 B + 0.96825 A 0 + 0 ,06052A S 

+ 0.003?8A 4 + 0.00024A 6 (79) 

The function B(y) of equation ( 74) is evaluated, at 
the nxrre -joaint s considered in table 6. 

The nine linear equations involved in the least- 
squares procedure are. to be obtained by evaluating (73) 
at the nine points y^ = k/8 , k = 0, 1, 2, ... 8.. 

Ifith the approximation of eq;uat ioii ( 17 the' matrix *f 
coefficients of the unknown parameters B and A 2n . iu 
this set of equations, corresponding. to. the " 1 eft— hand 
side of equation ( 21 ), is- obtained as the ' differ ence be- 
tween two matrices, the first of which is given directly 
by the left— hand side, of equation ( 19 : ) . In the second 
matrix, which represents the coef f ioient s ' of the values 
of the third term on the left-hand ' s ide of equation (73), 
the elements in the" k—th row are ' the products " of ' the 
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.coefficients of the parameters .. in equation (79) by the 
quantity 6 (y fc ). 

The coefficients of the required nine equations aro 
thus found, in the form 



B 


A 0 


A 2 


A 4 


A 6 




— s-3449 0 


.75651 


-.12069 


-.02512 


-.01182 






-.12560 


.76682 


-.09526 


-.02897 


-.01303 




i 


.04424 


.82172 


-.01896 


-. 03488 


-.01730 




i 


.19912 


' .83092 


.07988 


-. 03090 


-.02530 




i 


. 35128 


.81698 


. 20800 


. 00536 


-.02794 




i 


.48752 


.77392 


.35474 


. 09908 


.00242 




i 


.59186 


.69876 


.50264 


.27201 


.12231 




i 


.63452 


.57466 


.61642 


.52371 


.41573 




i 


.40518 


. .31474 


=.47045 


.6 329 2 


.76172 




i 



(80) 



The coefficients of the final set of linear equations 

are then found, by taking the dot product of the auxiliary 

matrix previously defined into both sides of equation ("80), 



B 


' A 0 




A s 


A 4 


A-6 - 






1.97915 


" 2.30991 


1 


.66370 


1.04882 


.74071 




3. 


45862 


2.30991 


6.55641 


1 


.75463 


.88374 


.54525 




8. 


71573 


1.66370 


1.75463 


1 


.45675 


1.00542 


.74857 




2. 


77811 


1.04882 


.88374 


1 


. 00542 


.84761V 


.71317 




1. 


67224 


.74070 


.54525 




.74857 


.71317 


.65352 




1. 


21165 



(81) 



and tho solution of this set of equations is determined, 
B =-1.45207 A 0 = 1.29708 A s ..= 2.59819' 

A 4 <= -2.38087 A e> = 2.03971 

In table 6 and figure 6 of the auxiliary lift function 
F(yJ - 



(82) 



rac 



it 



i l , given the. expression ■ 

'(7) - i'i(y)-. ~ — P(y) 

rr n. — '2 Ya log a 



(83) 



and the function 



J(y) 



cj are presented numerically 



o*(y) m 

and graphically. For the purpose of-compar ison there is 
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included in figure 6 a curve representing the variation 
6f the function — c 7 for the sane wing without cut-e-\ou£, 

EX » 

as determined in a preceding section. 



COU0LUSIO1T 



The explicit procedure given in this paper, using a 
five— tern approximation and a nine— point weighting system, 
should give results with a degree of accuracy comparable 
to the accuracy of measurement of the physical data of 
the wing in all practical cases. More nearly exact solu- 
tions could be obtained by retaining a greater number of 
terms or by using a more elaborate weighting system. 

While the present method is not a method of successive 
approximations, an -indicat ion of the accuracy attained can 
be found in the following two ways J 

1, As mentioned in the text, if the computed values 
of the parameters specifying the lift function P(y) 

are substituted i-n. the left— hand sides of the nine linear 
equations represented by formal equation (2l) the re- 
sultant values of the angle of attack corresponding to 
the function 3? can be compared with the corresponding 
prescribed values of the angle of attack. 

2. It can be shown that if the last term in the ap- 
proximation to the function S"(y) is not retained, the 
final set of four linear equations given by the least- 
squares procedure for the determination of the four re- 
maining parameters is obtained from the corresponding 
set of five equations in the original five parameters by 
deleting the coefficients of Ag and omitting the 'last 
equation. If the resultant set of equations is solved, 

a comparison of the lift distributions corresponding to 
the four— term and five— term approximations will afford 
an estimate of the effect of retaining additional terms. 

The complete solution of a problem, including the 
tabulation of values of the functions I*(y) and a(y), 
can be carried out in less than two hours by a competent 
cqmput ing-^machine operator. 

Massachusetts Institute of Technology, 
Cambridge, Mass., February 19 43. 
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TABLE 1.- LUI DISTRIBUTIONS FOB,' RECTANGULAR WING- ( — = m 



if 




1 
ma 


m » - . 







Present .Solution. 


Glauer \ 


*olut ion 


ft, — dim 


a .— aCmf y I 


a = a^y 2 . 


a = cticjt. 


a = arp[y] 


0 


0.8288 


■ 0.1673 


0.0666 


0.832 


0.132 . 


. 125 


.8262 


. .1898 


. 0757 


.830 


.157 


. 250 


.8243 


. . 2462 


.1075 


.822" 


.222 


.375 


.8169 


■..3171 


.1562 


.8 09 


,307 


.500 


.7932 


. .3820 


.2121 


.786 


.385 


.625 


.7478 


.4272 


Z2686 


.752 " 


.435 


.750 


.6828 


.449 8 


. 3236 


.695 


.451 


,875 


.5826 


.4362 


.3567 


.577 


.421 


.9 325 


.4736 


.3760 


.3235 


.450 


.355 


1 


0 


0 


0 


0 


0 


I ABLE 2.- 


- COMPARISON 


OF ANGLE Or 


ATTAOK a 


CORRESPONDING- TO 



TEE LIPT DISTRIBUTIONS OP TABLE 1 ¥IIH PRESCRIBED' 



ANCLE OP ATTACK a 



y 


a 




. * 


a = 


■■ a%[ y '] 




a = a T y 3 


a/ 










a/dfcrji 


Present 
Solution 


Glauert 
solution 


Present 
solution 


Claaert 
.gdution 


Present 
solution 


0 


0.997 


■1.000* 


1.0 


0.054 


0. 000 


o.oov 


- -0.002 


0 


.125 


.986 


•1 . 000 


1.0 


.109 


.050 


.125 


. 004 


.016 


. 250 


1.003 


1 . 000 


i.-o 


. 233 


.19 6 


.250 


.065 


.06.2 


. 375 


1.022 


1.003 


1.0 


.385 


.376 


.375 


.159 


.141 


,500 


1.018 


.995 


1.0 


,521 


.544 


.500 


. 265 


.250 


.625 


.982 


.995 


1.0 


.620 


.646 


.625 


.376 


.391 


.750 


.954 


1.005 


1.0 


.711 


*729 


.750 


.526 


.562 


.875 


1.030 


1.020 


1.0 


.894 


»851 


,875 


.79 0 


.766 


1 . 


.952 


.692 


1.0 


.974 


.861 


1. 


.962 


1. 



NAG A" Technical ''Note No. 925 



.SABLE 3.— i'AUXIJilAlty DATA .FOE TAPERED WING- 



3? 


c*(y) 


a(y) 




a*(y) 


0 


1. 


1.0 


2.17796 


1 .85269 


.125 


.9375 


1.0 


2.13103 


1.94784 


.250 


,8750 


1.0 


1.98288 


1.94088 


.375 


.8125 


■ 1.0 


1.69962 


. 1.76657 


.BOO- 


.7500 


• 0 


1; 06 538 


1.09523 


KS' 


.6875 


0 


-.5469 2 


.47025 


.750 


.6250 


0 


.26128 


. 09 278 


.875 


.5625 


0 


-.08 38 0 


-.17629 


1 . 


0 ■ 


0 


o- 


-.32527 



TABLE 4.- LIET D ISTEIBUTION JOE TAPERED VINO 



Present 


s olut ion 


Pear s on 


solut ion 


¥ . 


c 


y 


c 


m Cr C l 


mc E Cl 


0. 


0.7729 


0, 


0.7576 


.1250 


.7077 


.1564 


.7033 


.2500 


.6588 


.309 0 


.6535 


. 3750 


.6162 


.4540 


.4451 


.4375 


.5463 


.4890 


.3589 


.5000 


' .3057 


.5225 


.2743 


.56 25 


.1619 


.58 78 


.1460 


.6250 


. .1101 


.7071 


.066* 


. .7500 


.0744 


.8090 


.0306 


.8750 


.0425 


-.8913) 


.•0328. . 


1. 


0 


.9511 


.0224 
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» 

TABJLE 5.- COMPARISON OS 1 ANGLE 02* ATTACK a CORRESPOND ISO- 
TO THE LIPT DISTRIBUTIONS OP TABLE 4 

i 

WITH THE PRESCRIBED ANGLE OP ATTACK a 



Present solution 


Pearson solution 


y 


3<7> 


a(y>' 


y 


co( y) 


a(y) , 


0 

.125 
.250 
.375 
.500 
. 625 
.750 
.6-75 
1.000 


1.157 
.936 
1.008 
1.100 
-.023 
-.058 
.029 
. 020 
-.036 


1.0 
1.0 
1.0 
1.0 

0 
0 
0 
0 
0 


0 

.1564 
. 3090 
.4540 
.4890 
.5225 
.5878 
.7071 
.8090 
.8910 
.9511 


1.043 
.9 47 

1.058 
.724 
.505 
.288 

-.010 
. 063 
.021 
.046 
,044 


1.0 
1.0 
1.0 

1-rO 

0 
0 
0 
0 
0 
0 















TABLE 6.- NUMERICAL DATA POR WING VITH CUT-^JUT 



y 


P(y) 


P(y) 


c 

mc R * 


1 

m" °l 


0 

.125 


1.70627 
1.5991 1 


0.44512 
.42638 


0.9339 
.9 234 


0.9 339 
.9234 


"."25 0 


1113364 


.3449 8 


;9?60 


T.9760 
\.7320 


.375 
.500 
.625 
,750 
,875 

1 


.62280 
.37225 
.20701 
.0949 7 
.0249 2 

0 


.05355 
. 02069 
-.00098 
-.01568 
-.02486 
-.02813 


1.0240 
1.0272 
.9781 
.8965 
.7659 

0 


.7680 
.7704 
.7336 
.6724 
.5744 

0 
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Figs. 1, 
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rig. 8b ,o 



.8 
.6 

.2 
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Pr«»»nt idllillon 
— Glauirt Solution 
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(b) Symmetrically linear angle of attack (a - a<i |y|). 
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(c) Quadratic angle of attack (o£ » o^y 2 ) . 
Figure 3.- (Concluded) . 
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Figs. 3,4 




4 - Lift distribution for a tapered wing with rounded tips and partial-span 
flaps (b2/S = 10). 
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Figs. 5,6 
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• ••• Pearson solution 
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Figure 5.- Comparison of angle of attack corresponding to lift distributions of 
figure 4 with prescribed angle of attack. 
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Figure 6.- Lift distribution for reotangular wing with cutout (b/cmax = m, b 2 /S 



